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Abstract 

A simple transformation converts a solution of a partial differential equation with a Dirichlet 
boundary condition to a function satisfying a Robin (generalized Neumann) condition. In the 
simplest cases this observation enables the exact construction of the Green functions for the wave, 
heat, and Schrodinger problems with a Robin boundary condition. The resulting physical picture is 
that the field can exchange energy with the boundary, and a delayed reflection from the boundary 
results. In more general situations the method allows at least approximate and local construction 
of the appropriate reflected solutions, and hence a "classical path" analysis of the Green functions 
and the associated spectral information. By this method we solve the wave equation on an interval 
with one Robin and one Dirichlet endpoint, and thence derive several variants of a Gutzwiller-type 
expansion for the density of eigenvalues. The variants are consistent except for an interesting 
subtlety of distributional convergence that affects only the neighborhood of zero in the frequency 
variable. 



* Current affiliation: Systat Software, Inc., Richmond, California 

^Electronic address: fulling@math.tamu.edu| URL: Ihttp : / / www . math . tamu . edu/~f ullin| 



I. INTRODUCTION 



Here we develop a technique for constructing a solution to a differential equation with a 
Robin (generalized Neumann) boundary condition when a solution to the same or a related 
equation with the Dirichlet boundary condition is available. The idea is surely not new. (In 
fact, our key formula (j3.18|) for the heat equation was published in 1891 0], and Tikhonov 
and Samarsky indicate how to solve the wave equation on the half-line by a method 
related to ours. In both cases the initial heuristic motivation is somewhat different from 
ours, but the resulting calculus is the same.) But we believe that it has not heretofore been 
systematically developed and exploited. 

We are primarily interested in the integral kernels (Green functions) that solve the wave, 
heat, Schrodinger, . . . equations associated with a given self-adjoint (usually positive) dif- 
ferential operator, say H, in the spatial variables. These functions are useful not only to 
solve the partial differential equation concerned, but also to obtain information about the 
spectrum and eigenfunctions of H. In particular, semiclassical approximation (or classical- 
quantum duality |9|) relates the eigenvalues of H (collectively) to the periodic orbits of the 
classical system whose quantum Hamiltonian is H (or of the geometrical optics of the wave 
equation of H). 

In the simplest cases (the Laplacian operator with zero- dimensional or at least flat bound- 
aries) exact solutions of the Robin problems can be found. These kernels manifest a certain 
nonlocal behavior in the time variable, which is reflected in the more familiar eigenfunctions 
and frequency-domain Green functions by a nonpolynomial dependence of the reflection co- 
efficients on frequency {u in (j3.15|) . for example). This property makes it difficult to solve 
the time-dependent problems by simple matching at the boundary — hence the utility of 
the Dirichlet-to-Robin transform technique. The phenomenon of time delay at the bound- 
ary can be understood physically by observing [i^ that in a wave equation the Robin 
condition models an elastic support at the boundary; in other words, in one dimension a 
vibrating string is attached at one endpoint to a discrete, massless spring with Hooke con- 
stant K. The string can exchange energy with the spring; this explains the surface energy 
and action associated with Robin boundaries in quantum field theory [13, EH IH; 
also the possibihty of delayed reflection of a physical impulse. 

After setting up the general formalism of the method in Sec. IHl in Sec. Illll we construct 
the Green functions for a variety of time-dependent problems on the half-line. For the heat 
and Schrodinger equations the extension to a higher- dimensional flat boundary is easy, and 
the analysis leads to a (possibly new) determination of the heat-kernel coefficient (e.g., |5|]) 
associated with /t" for any n. Sec. II VI discusses the hopes for extending the method to less 
elementary models, where only approximate solutions can be expected. Sec. |V| treats the 
wave equation on an interval with one Robin endpoint; therefrom, the eigenvalues and the 
local spectral density are recovered from sums over the periodic and closed orbits of the 
problem. In fact, we stress that there are quite a few different ways of arriving at such 
sums and it is not always obvious that the results are the same. Most notably, because 
the series are not absolutely convergent, their behavior is sensitive to reordering of the 
terms; we show, however, that this problem is significant in practice only at very small 
frequencies and amounts in principle to a delta-function ambiguity at zero frequency (see 
Appendix^. This observation probably has implications for more general problems. Apart 
from that phenomenon, we demonstrate agreement among several variants of the periodic- 
orbit expansion. A later paper 0| will derive the Casimir energy for a Robin plate [2^ by 
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the frequency-cutoff method (cf. Sec. IIII El and [3, 



Notational remarks: 

1. For a half-space {x > 0} we write a Robin boundary condition at x = as 

^(0) = +KU{0). (1.1) 

Note that ^ is the inward normal derivative. Thus k equals —7 in the notation of 
jisl. in the notation of [2^, —S in the notation of and +n in the notation 

of 0, Q, . With this sign convention, k > is the physically more plausible case. 
(In a heat conduction problem, for example, k < represents heat flow from a colder 
medium to a hotter one.) When k > all the operators studied in the present paper 
have nonnegative spectrum. 



2. Like j23[ we use to stand for "irrelevant transverse dimensions" although the no- 
tation X|| (for "dimensions parallel to the boundary") would be equally logical. 

3. In this paper we find it convenient to use G as a generic notation for all Green functions, 
rather than introduce separate letters for heat kernels, wave propagators, etc. The 
meaning of G is stable within each subsection. 

4. When a Dirac delta function appears at an endpoint of an interval of integration, 
its interpretation is ambiguous by a factor |. We adopt the convention that the 
density of (nonnegative) eigenvalues, Ek = Uj?, with respect to integration over the 
eigenfrequency lu = \/E, is 

00 

p{uj) = J2^i^-^k) (1.2) 

k=l 

even when is an eigenvalue. The reader may note (especially in Sec. IV B|) some 
factors I that are not fully explained but are clearly necessary to produce the right 
answer. A rigorous treatment of such issues may appear elsewhere |l2| . 

5. All equations are in natural units, where time has units of [length]^ in heat and quan- 
tum problems but of [length] in wave and vacuum-energy problems. The only "cou- 
pling constant" that appears is k, with units [length]"^. The only dimensionless small 
parameters are ratios of the lengths k"^, t (or u)~^), and (in Sec. |^ L. Introduction 
of h into this context would merely be obfuscatory. 

6. We count dimensions nonrelativistically (e.g., a vibrating string obeys the one-dimen- 
sional wave equation). 



II. A DIRICHLET-TO-ROBIN MAPPING 
A. The Robin-to-Dirichlet map, T 

Consider functions /(x, x_l) defined on the half-space R^J_ = {(x, x_l) G R'^: a; > 0, x_i_ e 
R'^^^} and satisfying reasonable technical conditions (made more precise below). We define 
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(for a given constant k G R) 

^/-l-/- (2-1) 

Our first lemma is a tautology, but will prove to be powerful: 

Lemma 1 / satisfies the Robin boundary condition at x = if and only ifTf satisfies 
the Dirichlet condition, T/(0,x_|_) = 0. 

Now allow / to depend on an additional variable, t, and consider differential equations 
of the general form 

VV - + Vi/ = 5/, (2.2) 

where the operator S does not involve {x, x^) or derivatives with respect to them; for exam- 
ple, 6 = d/dt gives the heat equation, and 5 = gives the d-dimensional Laplace's equation. 
Let S stand for the natural domain of each problem: R^J_ for Laplace's equation or the 
eigenvalue problem (Sec. I11D|) . R'^"^^ for the wave equation, and R++^ in problems where t 
is inherently positive f Sees. IXTm and ITll E|) . 

Lemma 2 If f solves V^/ = 6f in S, then so does Tf. 

Lemma |21 follows immediately from the commutativity of T with and 6. It therefore 
hinges on the facts that has constant coefficients, k is a constant function, and the 
boundary is flat (so that the normal differentiation is d/dx). Nevertheless, we expect our 
ultimate construction to be useful in more general problems, as explained in Sec. IIVI 



B. Construction of T ^ for k > 

Given a function g, we wish to construct an / such that g = Tf. To make / unique we 
must impose a supplementary condition. 

The differential equation to be solved (with x^ suppressed for notational simplicity) is 
f'{x) — k/(x) = g{x), whose general solution is 



fix) = e^^ 



/ e-^'g{s)ds + Ce^''. (2.3) 
Jo 



It is natural to choose the solution with minimal growth as x oo. li g obeys a reasonable 
growth condition (such as boundedness) , then one can set 

POO 

C = - e'^'g{s)ds (2.4) 
Jo 

to cancel the exponentials and get a solution of similarly reasonable growth: 

POO 

fix, xx) = T-'g = -e^- / e-'-'gis, x^) ds 

J X 

e-^'gie + x,x^)de. (2.5) 



That is, if the domain of T is suitably restricted, an inverse operator exists and is given by 
formula ()2.5|) . A short calculation verifies that commutes with V^, and so we have a 
converse to Lemma El 
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Lemma 3 If g solves V^/ = Sf in S, then so does T ^g. 

The condition (yf(0,x^) = is not used in the proofs of Lemmas |21 and 01 But now put 
them together with Lemma ^ 

Theorem 1 If g solves the Dirichlet problem for V^/ = 6f in S, then T^^g solves the 
corresponding Robin problem (with the given k), and vice versa. 



C. Construction of T for k < 

This time it is convenient to treat the case n = 1 thoroughly before introducing the com- 
phcation of transverse dimensions. In that case y^{x) = e**^ is a normahzable eigenfunction 
of the problem; it satisfies Ty^: = 0, and therefore T can't be invertible. The most convenient 
growth condition is to require the functions to be square-integrable, so that a generalized 
inverse can be defined in the Hilbert-space orthogonal complement of . 

Suppose, for example, that u{t,x) is to satisfy the heat equation with initial data f{x). 
The part of the solution proportional to y^, can be written down immediately as e'^ ^Pf 
where 

^/ = 7rn§^*(^) (2.6) 
\\yA\ 

is the orthogonal projection onto y^ . Then the full solution is u{t,x) = u±{t,x) + e'^^^Pf, 
where u±{0,x) = (1 — P)f and u±(t,x) remains orthogonal to y^: at all t. Our task is just 
to construct u± . 

Accordingly, we now treat T as an operator in the Hilbert space range(l — P) and demand 
that the solution (j2.3|) lie in this space — i.e., {y*\f) = 0. It is easy to see that then (j2.4p 
is replaced by 

POO 

C = - e^'^'g{s)ds = -{y.\g), (2.7) 



and 



f{x) ^ T-'g = e'^^ 



e-^'g{s)ds- / e+"^^(s)rfs 
L^o Jo 



(2.8) 



Remarks: (1) In both cases one can write C = — e~^'^^'^ g (s) ds . (2) Because of the 
exponential decay of y^ , ()2.8|) makes sense for many functions that are not square-integrable. 
Furthermore, both ()2.5|) and ()2.8p can be applied to certain distributions by duality, and 
that will be done without comment in later sections. 

When n > 1 the kernel of T, as an operator in L^(R|j|_), consists of products of e^^ with 
square-integrable functions of x^. It is still true that the appropriate formula is 



f{x,^^)^T-'g = e^^ 



'^''g{s,:s.±) ds - / e'^'^''g{s,:s.±) ds 



(2.9) 



and that the domain of T is characterized by the requirement that e'^^f{x, x^) dx = (for 
all xx). This conclusion can be justified by performing a Fourier transformation in x^ and 
applying the reasoning above to each Fourier component separately, or merely by verifying 
that (j2.9j) satisfies all the necessary conditions. Thus Theorem HJ suitably interpreted, 
applies to both positive and negative k. 
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D. The eigenfunctions 



As a first "application" of Tfieorem ^ we clieck tliat it yields tfie correct eigenfunctions 
(solutions of V^/ = —u'^f) of the Robin problem. Ignoring the inert transverse dimensions, 
we have the Dirichlet eigenfunction g^j^x) = sin(ci;x) for each u > 0. Applying (j2.5|) or (j2.9|) . 
as appropriate, gives 



-UJ 



uj' + k' 



cos[uJX) H sm[ujx) 

UJ 



(2.10) 



which agrees up to normalization with the Robin eigenfunction given in standard references 
0, H^. (Of course, this is the hard way to reach an elementary result.) In Sec. IIIIBl we shall 
obtain the normalization by our methods as well. 



III. APPLICATIONS TO ELEMENTARY GREEN FUNCTIONS 
A. Generalities 

Henceforth we restrict attention to k > 0. 

A Green function (integral kernel) associated with a Dirichlet problem in S typically has 
the "image charge" form 

Goix, x±, y, y±) = G{x -y,x±- yx) - G{x + y,x±- y±), (3.1) 

where G{x — x_|_ — y_|_) (which is even under interchange of x and y as well as translation- 
invariant) is the corresponding Green function for all of R*^ (and the time variable, if any, 
is momentarily suppressed in the notation). 

These Green functions represent operators that are functions of A and hence commute 
with T. Therefore, in operator notation, Gr = T^^GdT should be the corresponding 
operator for the Robin problem. It is understood that the action of a Green function on a 
function is 

G/(x,x_l)=/ dy I (iy_LG(x,x_L,?/,y±)/(y,y±). 

Therefore, the multiplication by T on the right is represented by the transpose (real adjoint) 
of T acting on the variable y: 

Tt = - — 
^ dy 

The multiplication by on the left is represented by the corresponding integral operator 
applied to the x variable of G. 

On a function of the form G{x — y), is equivalent to . On a function of the form 
G{x + y), passing from y to x leaves T'^ as T\ which can also be written as 

= -T - 2k. (3.2) 

Since T commutes with G, the effect of the similarity transformation is to leave the first 
(direct) term of unchanged, while the second (reflected) term changes in a rather simple 
way, stated in a mixed operator/function notation (with transverse variables suppressed) in 
the next lemma. 
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Lemma 4 In the context of iy.l]} . 

(a) T-'G{x-y)T = G{x-y), 

(b) -T-^G{x + y)T = +G{x + y) + 2kT~^G{x + y) . 

Putting the two parts of the lemma together, we get our principal working equations for 
the rest of this section, (j3.4|) - (j3.5|) . 

Theorem 2 When a Dirichlet Green function has the form 0^.1}) . the Green function for 
the corresponding Neumann problem is 

Gn{x, x±,y, y±) = G{x -y,x_L- yx) + G{x + y,x±- y±), (3.3) 

and the one for the Robin problem is 

Gnix,x^,y,y^) = GN + A^G (3.4) 

where 

A^G{x, x_L, y, yx) = 2kT^^G{x + y, x_l - y±), (3.5) 

B. The wave equation in one space dimension 

Consider the wave problem 

The well known d'Alembert solution corresponds to the Green function 

G{t, x-y) = \[6{x -y-t) + 6{x-y + t)]. (3.7) 

From (|3.ip . therefore, 

Goit, x-y) = \[5{x -y-t)+5{x-y + t)]- \[5{x + y - t) + 5{x + y + t)], (3.8) 

where the last of the four terms is relevant only for t < since x and y are positive in the 
physical region. Thus, by Lemma IH we are interested in 

- \ T-^6{x + y-t)T = ^6{x + y-t) + kT-^6{x + y-t), 

POO 

T-^6{x + y-t) = -e'^^ / e-^'6{s + y - t) ds 

J X 

= -e"("+^-*)0(t-y-x), 
where 6 is the unit step function. So, finally, for t > and k > we have 
Gji{t, X, y) = l[6{x -y-t) + 6{x-y + t) + 6{x + y-t)]- /te"(^'+^"*)^(t -x-y). (3.9) 
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y X y X 



FIG. 1: (a) For the wave equation with a Robin boundary, there is exactly one delayed path 
(solid curve) from space-time point {0,y) to a generic point (t, x). With a Neumann or Dirichlet 
boundary, no delay occurs and generically the reflected path (dashed curve) misses {t,x). (b) For 
the Schrodinger equation, there is exactly one path without delay from (0, y) to {t, x) (dashed 
curve); it has constant speed vq . When delayed reflection is possible, there are inflnitely many 
delayed paths (solid curve) with all speeds v > vq . 



That is, the solution of the wave equation with Robin boundary condition and initial 
data u(0, x) = /(x), f{y) being interpreted as when ?/ < 0, is (for t > 0) 

n(t, x) = i[/(x - t) + fix + t) + fit - x)] - KOit -x) f ^ e-^^'-^-y^fiy) dy. (3.10) 

As usual in one-dimensional wave problems, the solution contains a right-moving unreflected 
pulse, fix — t), and a left-moving pulse, fix + t), which reflects off the boundary when t = x 
as a new right-moving pulse, fit — x). In addition, in the Robin problem there is a smeared- 
out term, which is perhaps more illuminatingly rewritten 

pt—x 

-neit-x) I e-^Jit-x-e)de. (3.11) 



Physically, this formula indicates that not all of the pulse is reflected immediately — the 
wave continues to leak out of the boundary with an amplitude that decays exponentially in 
the time delay e. There is not only an "echo" but also "ringing". Otherwise put: In the 
Dirichlet and Neumann problems, the signal at x caused by a source at ?/ is a sum over all 
the unit-speed "classical paths" from (0,?/) to (t,x), including possible reflections; there is 
at most one such path (for most points none). But in the Robin problem one must also 
integrate over all the paths that travel from iO,y) to the boundary with unit speed, stay 
there for a time e, and then go with unit speed to (t, x); for t > x + y there is always exactly 
one such path (Fig. HJ^a)). Interference between the two terms of (|3.1U|) is responsible for 
conservation of energy. 

Formula ()3.9p is the central result of this section. It can easily be generalized to negative 
t and also to problems with nonzero initial data for For later reference we note that 

/*00 /"OO 

/ A^Git,x,x)dx = -K e'^(^^-*)e(t-2x)rfa; = i(e-"*-l). (3.12) 
Jo Jo 

Finally, we can recover the eigenfunctions, complete with normalization (relative to the 
Lebesgue measure duj). The generalized eigenfunction expansion 

POO 

Gnit,x,y)= / cosiujt)i!^ix)i!*M duJ (3.13) 
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must hold, hence 

2 



7r 



^oj{^WM = - GR{t,x,y)cos{iut)dt. (3.14) 



From ()3.9|1 one gets for the right-hand side of ()3.14j) 

— ^-^ — - [a;^ cos (jjx cos ujy + fi;a;(sin uox cos ouy + cos sin ujy) + sin ujx sin a;?/] , 



which can be factored as '^/'^(x)-?/'* (y) with 

i'ioix) = {uj'^ + k'^)"^^'^{uj cos UJX + K sin ujx) (3.15) 
2 

— sin(co'x + 0), 

where uj uj k 

(j) = ta.n~^ — ; sm(j)=—^===, cos (p = —^=== . (3.16) 



1^ VuF+i^ \fuF^ 

This normahzation agrees with that in the treatises 0,113]. 



C. The heat equation 



The Green functions for the heat and the Schrodinger equation are essentially the same 
algebraically, one being an analytic continuation of the other. In this subsection we derive 
the heat kernel by the Dirichlet-to-Robin transformation and compare with the result of a 
direct eigenfunction expansion. In the next subsection we recast the Schrodinger kernel as 
a "sum over classical paths". 



The heat kernel for R is 



G(t,x,x^,y,y^) 



{A-ntyi'- 



exp 



[x - yf + (xx - y±) 
4t 



(3.17) 



The heat kernels for the Dirichlet and Neumann problems in the (i-dimensional half-space 
are then given by ()3.1|) and ()3.3p . These functions are products of the corresponding one- 
dimensional functions of x and y by the free id — l)-dimensional heat kernel, which is 
unaffected by the operations T and (for constant k). Therefore, in what follows we look 
only at the one-dimensional kernel to streamline the notation. 

According to Theorem |2l the Robin heat kernel equals the Neumann kernel, Gj^ , plus 



A^G = iKT-^Git, x + y) = -2k- 



(47rt)i/2 

This expression can be rearranged into 

[s + y + 2Ktf 



/■oo 


{s + yf 


/ e"""" exp 

J X 


At 



ds. 



(3.18) 



K 



{TTty/- 



exp 



At 



ds = -fi:e'^(^+2')e«'*erfc ( + K^ft] . 

\VTt J 
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Here erfc is the complementary error function, denoted by 1 — $ in [18[. Thus, finally, we 
have (for k > 0) 



GR{t,x,y) 



(47rt)V2 



|exp 


{x-yf 
At 


+ exp 


{x + yf 
At 


} 










erfc I 



Formulas ()3.18|1 and ()3.19|1 (which are not new [^, 0|) are the two key results of this 
subsection; (|3.18p has a path-sum interpretation, which is best postponed to the next sub- 
section. 

From ()3.15j) . the eigenfunction expansion of the Robin heat kernel is 



GR{t,x,y) = - 
vr 



cos(co'x) H smiujx) cosiuy) H sm(uy) 



(3.20) 

The integrals can be evaluated by formulas (3.954) of [18] (printed incorrectly in some earlier 
editions); eventually the same result (j3.19|) is obtained, but the calculation via is quicker. 
An integration by parts shows that 



^^Git, X, x^, X, xj_) dx = (Ant) 



(e" *erfc(/s:v^) - l). 



(3.21) 



Expanding ()3.21|1 as a power series in k or t^/^ (see (8.253.1)]), 

^^^)^_J_^2^{KVtf^+^ 



i(e""erfc(KVi)-l) = iE^-4=E 



K\/t 



(2j + l)!! 



2 



+ 



+ 



^ , , , (3.22) 

37? 4 ' ^ ' 

reproduces the known contributions jl, HH of a (flat, constant-^) Robin boundary to the 
usual heat-kernel trace expansion and extends that information explicitly to all orders. 
Perhaps more interesting is that (j3.2H) is a rare example of a heat trace known exactly 
for all t, thereby leading to spectral densities known exactly for all u;; we return to this 
point in Sec. IV El 



D. The SchrSdinger equation 



The Green function for the time-dependent Schrodinger equation, = — V^m (in units 
where h = 1 and m = |), also known as the quantum propagator, is obtained formally by 
replacing t by it in all the equations of the previous section. It can be seen that this takes the 
variable to the boundary of the domain of analyticity where all the integrals are meaningful. 
More interesting than the resulting formula in terms of the analj^ic continuation of erfc is 
the counterpart of the prior formula ()3.18|) . 



GR{t,x,y) 



{Amty/^ 



exp 



■ {x - yf 
At 



2k 



+ exp 
1 



.{x + yf 



At 



(47rit)V2 



poo 


'.{s + yy 


/ e"""" exp 

' X 


r 4t J 



ds, (3.23) 
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which admits an interpretation as a sum over paths. In analogy with the previous discussion 
of the wave equation, one would expect a source at (0,y) to influence the solution at (t,x) 
along the direct path between the points (the contribution of the free quantum kernel), and 
along the path that bounces elastically off the boundary (the image term of the Neumann 
solution), and possibly along paths that hit the boundary, stay there awhile, and then return 
with the same energy; the problem is to show that the integral term (j3.18p lends itself to this 
last interpretation. The difference from the wave case is that in nonrelativistic mechanics 
the paths may have any speed. Therefore, for any two points there always exist a direct 
path and an echo path, and also infinitely many paths of the delay type (Fig^b)). 

Recall first that for a freely moving particle the action of a trajectory segment of length 
L, speed v, and time t is 

^ vL vH 

S = — = — = — . 3.24 

4t 4 4 ^ ^ 

The action functional is additive over segments, as is clear from the last form given, in 
which f^/4 is the constant (kinetic) energy of the orbit. So the quantities {x =F y)^/4 that 
appear in the exponents of the direct and echo terms of (j3.23|) are the total actions of the 
corresponding paths. 

In the integral term in (|3.23|) we make the usual change of integration variable s = x + e 
and the further substitution 

e = vt-{x + y) (3.25) 

(thereby defining v) and also define 

_ x + y _ x + y 

Vq = , u = t . (3.2o) 

t V 

The term becomes 

A^G = . / e-"*(''-"°)e™'*/^ dv (3.27a) 



(47r2t)V2 

2Ktl/2 



oo 



(47ri)i/ 



e 



-Kuv ^i[u(x+y)+v'^u]/A (3.27b) 



(where u depends on v). This equation has the following physical interpretation: v is the 
speed of the "particle" as it travels from y to the boundary and again from the boundary 
to X. Thus v{x + y) /A is the action of those two trajectory segments, {x + y)/v is the 
time consumed by them, and therefore u is the remaining time, which the particle spends 
somehow attached to the wall. The term v'^u/4, the product of this time with the energy of 
the orbit, is the action associated with this sojourn at the wall. Each such orbit contributes 
to the propagator with an amphtude — 2K)f:e~'^™ times the usual amplitude, (47rit)~^/^. Also, 
e/t = V — Vq is the difference (necessarily positive) between the speed of this orbit and that 
of the echo orbit, a limiting case. 

Remark: The discovery of the action expression A = w^m/4 (and hence the correct relation 
between e and v) was guided by the principle that the total action should be v{x + y)/4: + A, 
where A depends only on the local physics at the boundary; that is, A could be a function 
of V, u, and (in principle) n, but must not depend in any other way on t, x, and y. That A 
turns out to be precisely the total energy times the elapsed time u was an unforeseen bonus, 
as was the fact that the exponential factor in the amplitude is likewise determined by the 
local physics. 
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Since the integral ()3.18p for the heat equation has the same structure, conceptually it 
also can be given a path-sum interpretation. Because of the diffusive nature of solutions of 
the heat equation, however, any trace of "classical" behavior is difficult to discern, even for 
sharply peaked initial data. 



E. Laplace's equation in one higher dimension 

The Green function that solves V^m + |^ = (with u bounded as t ^ +00) for given 
boundary data u{0, x, x^) on the hypersurface t = is called the Poisson kernel or cylinder 
kernel for the spatial geometry concerned. Its limiting behavior as t | can be used to 
determine the vacuum (Casimir) energy density of a scalar field in that geometry |l5l |. 
When d = 1 and the initial hypersurface is the whole line, the Poisson kernel is 



Git,x,y) = - 



TT + (x — yy 

[Higher dimensions will be treated elsewhere [l7j.) 
In this problem the analogue of ()3.18|) is 



ds. 



Let w = s + y and perform a partial-fraction decomposition: 

1 1 



AC. 



TT 



x+y 



w + it w — it 



dw. 



By [18, (3.352.2)] one gets (see ^M)) 

AG. = e'^(^-+2^) [e'^tEi{-K{x + y + it)) - e-'''' Ei{- k{x + y - it))] , 



which can also be written 

AG. = ^ e''(^+*') Im [e"*"* Ei(mt - k{x + y))] (for t real). 

To get the complete cylinder kernel for the Robin problem, add 

t t 



a 



N 



1 

7T 



+ 



t^ + {x-yY t^ + ix + yY 



(3.28) 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



The eigenfunction expansion analogous to ()3.2()j) just has e in place of e From 



there a lengthy calculation using [18|, (3.354.1,2) and (8.233.1)] verifies that Gr — Gm = AG. 
as given in ()3.30p . In summary, this problem is very similar to the heat problem in Sec. IIII C\ 
but with different special functions appearing. 

The cylinder kernels discussed here solve Laplace's equation with nonhomogeneous Dirich- 
let data at t = 0; they suffice for calculating total Casimir energy and for calculating vacuum 
energy density when the latter is defined with the value of the "conformal coupling constant" 
set to C = i (j4 ll^. Obtaining the energy density for other values of ^ requires the Green 
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function for the problem with Neumann "initial" data, §j(0, x,x^). For the whole real line, 
that kernel is 

(3.33) 



Again, Theorem |21 can be implemented exactly to solve the temporal Neumann problem with 
the Robin condition at the spatial boundary, x = 0; the relevant integral is [18|, (4.337.1)], 
and the result is 

—{\Yi[e + {x-yf]-\Yi[e + {x-yf]} + -e^^''+y^ Re \e-'^' 'E.iimt - k{x + y))] . (3.34) 

(Note that the first term is , the kernel for the Dirichlet homogeneous spatial boundary 
condition; the reflection term in has been overwhelmed by an identical term in A^G 
with a factor —2.) 

Remark: In keeping with the well known "Green's identity" structure of the solution 
formulas for second-order elliptic boundary-value problems, the cylinder kernel for the non- 
homogeneous temporal Dirichlet condition is the t derivative of the corresponding cylinder 
kernel for the temporal Neumann condition. This relationship is easily checked for both 
pairs, (jSSHll/dSSni) and (nOTll / (n01|l . 

Just as the Schrodinger kernel is an analytic continuation of the heat kernel, replacing t 
by it in a cylinder kernel leads to a certain fundamental solution of the wave equation. More 
precisely, in ()3.28|) and ()3.30p (where t was positive) one should replace t by i(t — iO), where 
now t can have either sign but the negative infinitesimal imaginary part is needed because 
a singularity is encountered on the real axis when \t\ > x + y. For the basic Green function 
in ()3.28|) the singularity is a pole, which has the well known decomposition 



t 



71 {x- yf 



]-[6{x - y -t) + 5{x - y + t)] + 

Z 2,71 



1 



X — y — t X — y + t 



. (3.35) 



The delta term is recognized as the d'Alembert Green function fl3.7|) . which solves the wave 
equation with given initial value and vanishing initial time derivative. The principal-value 
term appears because the asymptotic condition on the Poisson kernel at t = +oo has evolved 
into a positive- frequency condition on the solution of the wave equation, hence an initial value 
of % that is a certain nonlocal functional of the initial value of u. In the language of quantum 



is (proportional to) the time derivative of the Wightman function, and 



dt 

field theory, 

the d'Alembert term is proportional to the time derivative of the field commutator. (Starting 
from the other kind of cylinder kernel would avoid the time derivatives.) 

All these statements have analogues for the Green functions of the Robin problem. From 
fl3.30|) we get as the correction to the Neumann Wightman function 



AG, = — e"("+^)[e-"*Ei(fi:(t 

TT 

Here Ei is defined by [Ti, (8.211.1)]: 

Eif^) = 



X 



y))-e'''Ei{K{t + x + y))]. 



dt 



(3.36) 



(3.37) 



with a branch cut on the positive real axis (which comes into play when \t\ > x + y). 
According to HI, (8.240.3)], Ei(a; T iO) = VEi{x) ± ztt. (See also [H (3.352.5)] and the 
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footnote on p. 228 of (H.) Thus ()3.36|) consists of a principal- value term, plus a jump term 
that equals —Ke'^'^^^'^~^^6{t — x — y) in the case of positive t (for instance). This last is 
recognized as the A^G found for the classical wave propagator in ()3.4|) . 

As in the quantum problem one can in these problems introduce the variables 

X + y + e X + V 
v = , u = t —, 3.38) 

t V 

and interpret them as the speed of a reflected path from y to x and the time delay of the 
path at the reflection point. All that changes is the way in which the part of the integrand 
coming from the original Green function depends on f ; in the case (|3.29|) one has 

= - — r e-*(--o)_^ (3.39) 

7^ Jvo + 1 

in place of ()3.27aj) . As v varies from the kinematic minimum, vq = {x + y)/t, to infinity, u 
varies from to the kinematically allowed maximum, t. These reparametrizations are not 
useful calculationally in the simple problems treated in this paper, but they are likely to 
become important in arriving at a physically correct approximation ansatz in more compli- 
cated problems where the relation between time displacement (u) and space displacement 
(e) is nonlocal. The interaction between the field and the Robin boundary should take place 
at the boundary, not throughout a spatial layer of size e ~ 



IV. STRATEGIES FOR ADVANCED APPLICATIONS 

One would like to extend the Dirichlet-to-Robin technique to problems where the bound- 
ary is curved, k is not constant, or — 5 is replaced by a differential operator whose 
coefficients depend on x. There are two obstacles to be overcome. First, one must have a 
valid Dirichlet solution, v, or Neumann Green function, Gn , from which to start, and in 
general even those elementary boundary conditions cannot be solved exactly by the method 
of images. That problem leads into the general subject of semiclassical (or other) approxi- 
mations and will not be further considered here. Second, in general T will not commute with 
the differential operator, and hence u = T~^v will not satisfy the same partial differential 
equation as v. There are two strategies one might pursue to get around this problem. 

The first (which is not our favorite) is to find a differential equation to be satisfied by the 
Dirichlet function v that will cause the Robin function u to satisfy the correct differential 
equation. Note that there is some freedom in how to extend the definition (j2.ip of T 
to the interior region, since k is given only on the boundary. For example, consider the 
wave equation on the interval < x < 1 with a second Robin condition at the right end, 
— 1^(0, 1) = ^'^(0, 1) (with negative sign because the normal derivative now points in the 
opposite direction). If v{t,x) solves the (doubly) Dirichlet problem on the interval, then 
T~^v obeys the correct Robin condition at x = but not at x = 1, unless k' = —k. Not 
surprisingly, k' = —k is a condition for the eigenvalues of the Robin problem to be the same as 
those of the Dirichlet problem (see j23l . (4.4)]); it is clear that the construction Gr = T^^G^T 
is not possible for any T unless the problems are isospectral. It is possible, however, to choose 
a function hi{x) that smoothly interpolates between /€(0) = k and fi;(l) = —n' and to define 
T accordingly so that Lemma Q holds at both endpoints. The solution for now is more 
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complicated than ()2.5|) or ()2.8|) but still elementary. Finally, the equation that v = Tu must 
obey to cause u to obey the wave equation is 

(9^t> d'^v ^dn du cPu 

dt"^ dx^ dx dx dx'^ 

^dn ^ d / du ^ 2\ (4 1) 

dx dx \dx / ' 

an integrodifferential equation in general. With luck one can choose k{x) so that the coeffi- 
cient of u vanishes, and then one has a slightly modified wave equation for v (with a spatial 
operator whose Dirichlet realization is isospectral to our Robin problem). For some values 
of k{0) and k(1), however, the k,{x) satisfying the required conditions has a pole inside the 
interval (0, 1), so that this construction fails. 

In higher dimensions such difficulty appears to be generic. Defining T throughout a disk, 
for example, requires extrapolating the normal vector to the boundary smoothly throughout 
the interior. Such a vector field must vanish somewhere, and at that point the first-order 
partial differential equation to be solved to construct becomes singular. 

We conclude that although this approach may allow some special problems to be solved 
by tricks, it is not promising as a systematic method. 

The second strategy is to exploit the locality of the boundary interaction and the classical- 
path structure of the field dynamics. Our expectation is that a locally defined Dirichlet-to- 
Robin transform tells how any solution locally refiects off a Robin boundary, and that 
this information can be combined with standard semiclassical technology in the bulk (and for 
Dirichlet and Neumann reflections from curved boundaries) to construct global approximate 
solutions. In the next section we implement this approach for what is probably the simplest 
situation, two parallel flat boundaries with empty Euclidean space between them. 



V. AN INTERMEDIATE APPLICATION: WAVE EQUATION AND SPECTRUM 
ON AN INTERVAL 

A. The problem 

Consider flrst the one- dimensional wave equation on an interval with a Dirichlet boundary 
at the right end and a Robin boundary at the left: 

0^ = 0^ iO<^<L, 0<t<oo), (5.1a) 
u{0,x) = f{x), _(0,x) = 0, (5.1b) 

du 

— {t, 0) = Ku{t, 0) (k > 0), u{t, L) = 0. (5.1c) 
ox 

(There would be no difficulty in principle in handling a second Robin condition at x = L.) 
The makeup of the solutions from right- and left-moving pulses makes this model particularly 
easy and instructive. As shown in Sec. IIVI it is not possible to obtain a solution by applying 
once and for all to the solution of the pure Dirichlet problem (although the most severe 
complications mentioned in Sec. IIVI are not present here). Instead, the transform will be 
applied repeatedly as each pulse strikes the Robin boundary. 
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FIG. 2: Image sources for the wave equation on interval (0, L) with Neumann boundary at and 
Dirichlet boundary at L. Filled circles indicate positive terms, open circles indicate negative terms. 
The numerical labels are for comparison with Fig. EJ 



B. The Neumann— Dirichlet case 

To start, recall what happens when k = 0. The method of images associates the generic 
point y in the interval with an infinite set of image points in the unphysical regions, as 
shown in Fig. where the open circles indicate the images that are weighted negatively. By 
d'Alembert's formula ()3.7p . therefore, the wave Green function is 

^ oo 

G^(t, x,y) = - J2 (-1)" [^(t + x-y- 2nL) + 6{t-x + y- 2nL) 

n=—oo 

- 5{t + x + y -2{n + l)L) - 5{t - x - y - 2{n - 1)L)] . (5.2) 

With this choice of indexing, the pulses relevant for t > are those with n > in the first 
three terms and those with n > 1 in the fourth term, and when we go on to the Robin 
problem, n will be the number of times the operator needs to be applied. 

The trace of the wave kernel is the Fourier cosine transform (with respect to u) of the 
eigenvalue density. Here it is 



/ G'7v(t,x,x)rfx = V (-l)"[M(t-2nL) 
Jo 



-\o{0<-l + {n + 1)L < L) -^9{0 < I- {n-l)L < L)] , 

where d{a < t < b) = 9{t — a)9{b — t) is the characteristic function of interval {a,b). The 
inequality in the last term is equivalent to 2(n — 1)L < t < 2nL, while that in the next-to-last 
term is equivalent to 2nL < t < 2{n + 1)L. Reindexing then shows that these two terms 
cancel. (See also ()5.25|) . however.) Thus one has 



/ G'7v(t,a;,a;)c/a; = L ^ (-l)"5(t-2nL), 



(5.3) 



a sum over all the periodic orbits (with lengths 2\n\L). Its inverse Fourier cosine transform 



IS 



Pn{^) = z I dt cos(tu) I Giy(t,x,x) dx 
Jo 



TT 



or °° / 1 \ 

^(_l)Wl__5„Jcos(2nLcu). (5.4) 



TT 

n=0 
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This sum can be evaluated by the Poisson summation formula as 

J- oo J- oo 



L 

7T 



2n6{-2TTk + TT + 2Luj) 

=— oo 



k=l 



(since only > is meaningful in the cosine transform). That is, the eigenvalues (or, rather, 
their square roots, the eigenfrequencies) are . . . , as expected. 

Remark: The cancellation of the contributions from "bounce" orbits (closed but not 
periodic) is an artifact of the mixed Neumann-Dirichlet boundary conditions. In general 
those orbits produce the "surface area" term in the Weyl expansion of the eigenvalue density 
(cf. [26]). For example, in the pure Neumann case the (—1)" would be missing from the 
formulas above, so that 

/•L oo ^ oo 

/ G{t,x,x)dx = L V 6{t-2nL) + - V e{2nL < t < 2{n + 1)L) 
Jo 2 

" n=— oo n=— oo 

oo ^ 

= L Y ^(t- 2nL) + - . (5.6) 

2 

n=— oo 

In the pure Dirichlet case the | becomes — |. In the inverse cosine transform the ±| produces 
±^S{u!) (cf. Appendix . which combines with a +^6{uj) in the Poisson sum to denote the 
presence or absence, respectively, of the eigenvalue at = 0. 



C. The wave kernel 

Turn now to the Robin case. The wave kernel Gji will have the same basic pulse structure 
as Gat , but every time a pulse reflects from the left boundary it will acquire a time-delayed 
component; these effects cumulate as indicated in Fig. El At the first step we need to know 
how to produce the pulse numbered 4 in the figure by Robin reflection of pulse 1. Since 
pulse 1 by itself is not a Green function. Lemma |3] and Theorem |21 do not apply directly, 
but they carry over in essence: 

Lemma 5 Let u{t, x) be any solution of the one-dimensional wave equation. Then 
UR^t, x) = u{t, x) + u{t, —x) + 2kT^^ [n(t, —x)] 

pOO 

= u{t, x) + u{t, -x) -2k e''^^"'^n(t, -s) ds 

J X 

is a solution satisfying the Robin condition. 
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y 



FIG. 3: Structure of the solution (|5.7j) of the wave equation on interval (0, L) with Robin boundary 
at and Dirichlet boundary at L. 

Proof: u{t, —x) satisfies the wave equation by virtue of the latter's reflection symmetry. 
Then the third term is also a solution by Lemma 01 It remains to check the boundary 
condition: 



TuR(t, x) = Tu(t, x) + Tx[u(t, —x)] + 2Ku(t, — x) 
hence TuR{t, 0) = 0. 

Remarks: (1) The lemma is to be used in cases where u represents a wave impinging 
on the Robin boundary. Then the rest of ur represents a wave reflected in the opposite 
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direction. The reflection is "causal" in the sense that the incident wave is not modified 
by the construction. (2) This construction in its simple form may break down in more 
general problems (e.g., Schrodinger's equation with a potential), because (a) u{t,—x) may 
not satisfy the differential equation, or even be defined; (b) may not commute with the 
differential operator. See remarks in Sec. IIVI 

Pulse 1 is simply the term ^6(t + x — y) in Gjy . Apply Lemma to get pulse 4: 



-S(t — X — y) — 2k 



oo -| 

e<^-')5{t -s-y)ds = -S{t - x - y) - 2Ke-"(*~^-^)^(t -x-y). 



Continuing inductively, one builds up the entire wave kernel as presented in the following 
theorem. 

Theorem 3 For t > in the wave problem i5.1]} the Green function is 

oo 

Gnit, X, y) = [\5{t + x-y- 2nL) 

n=0 

- KLl_^{2K{t + x-y- 2nL))e~'^(*+^-^-2"^)^(t + x-y- 2nL) \ (5.7a) 

oo 

+ Y^{-lY\\5{t-x + y-2nL) 

n=0 

- KLl_^{2K{t -x + y- 2nL))e-^^'-''+y-^''^^e{t -x + y- 2nL)\ (5.7b) 

oo 

+ J](-l)"+' [|5(t + x + y-2{n + 1)L) 

n=0 



- kLI_^ {2K{t + x + y-2{n + 1)L)) 

^ ^^n{t+x+-y-2{n+l)L)Q^^ + 3; + ^ _ 2(n + l)L) 



(5.7c) 



^Y^^_l)n-ih5{t-x-y-2{n-l)L) 



n=l 



kLI^^ {2K{t -x~y-2{n - 1)1)) 



X e 



^^^e{t-x-y-2{n-l)L) 



(5.7d) 



Here L\_^ is the Laguerre polynomial fl3, (8.970.1) and (8.971.2)] 



X) 



-L'Jix 



{L\ ^ 0) 



(5.8) 



Remark: In terms of the labeling in Fig. ^ term ()5.7aj) comprises pulses 4n + 1, ()5.7b|) 
pulses 4n + 2, ()5.7c|) pulses 4n + 3, and ()5.7d|) pulses in. (Note that n starts from 1 in ()5.7d|l . 
but from in the other three terms.) 

The proof is a straightforward induction guided by Fig.|21 Robin reflection (by Lemma 
of pulses 4n + 1 and 4n + 3 produces pulses 4(n + 1) and 4(n + 1) + 2, respectively. Similarly, 
Dirichlet reflection at x = L of pulses 4n and 4n + 2 yields pulses 4n + 1 and 4n + 3. 
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D. The wave trace 



Let us now take the trace of ()5.7|) . The terms ()5.7a|) and ()5.7b|) give identical contribu- 
tions, totaUing 

/ Gn{t,x,x)perdx = N + P, (5.9) 
Jo 

where 

oo 

A^^L J](-l)"(5(t-2nL) (5.10) 

n=0 

coincides (for t > 0) with the result (j5.3|) found previously for the Neumann boundary 
condition and 

oo 

P = -2kL ^(-1)"L;,_i (2K(t - 2r2L))e~"(*-2"^)^(t - 2nL) (5.11) 

n=l 

is the additional contribution in the Robin problem of the periodic orbits. Note that both 
these terms carry an overall factor L — they came from an integrand independent of x 
— which is lacked by the remaining terms, which come from the bounce orbits and reflect 
the latter's close association with the boundaries rather than the global geometry. The 
contribution of the bounce orbits simplifies to 

L 

GR{t,x,x)i,oudx = A + B, (5.12) 



where 
and 



A = lie-^' - 1) (5.13) 



B = J2Y1 H — - 2nL)'"e-"(*-'"^)^(t - 2nL), (5.14) 



n=l m=l 



c{n,m) = ^(-l)«-i2^-i^!^('''j. (5.15) 

j=rri ^ \3 J 

Remarks: (1) For fixed t all sums encountered so far are finite, so there is no issue of 
convergence or term ordering. This property will be lost at the next step! (2) As expected, 
all terms except vanish as k ^ 0. As k ^ +cx3, it can be shown that B ^ Q, A ^ — |, 
and 

oo 

P^L5^(l-(-l)")5(t-2nL); 

n=l 

thus 

/ (jj^'(t, X, x) dx ^ L 5{t — 2nL) — as k — > oo, (5.16) 

which, as previously remarked, is the correct formula for the Dirichlet problem. (3) The 
critical dimensionless parameter of this model is Lk. So far all the calculations are exact, 
so it has not been necessary to assume Lk either large or small. (4) We omit the lengthy 
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calculation leading to (B), except to mention that it involves definite integrations over 
intervals of the form t — 2(n+l)L < s < t — 2nL, after which the index n in the contributions 
from lower limits of integration need to be shifted relative to those from upper limits to 
combine terms. 



E. The eigenvalue density 

We have already calculated pn ()5.4|) . the inverse cosine transform of A^. That of A is 

P«,avM = - - |5M («:^0). (5.17) 

(Contrary to appearance, this object does approach the Neumann limit of as k ^ 0, 
because the first term converges to ^S^u) in the distributional sense. The corresponding 
spectral-staircase formula, 

av(t<^) = - tan"^ - - ^ for > 0, (5.18) 

looks less anomalous.) Because A is recognized as the wave trace (j3.12p associated with 
the Robin boundary sitting in infinite space, one can interpret p^.av as the contribution of 
a Robin boundary (relative to the Neumann base case) to the "averaged" or "smoothed" 
spectral density in dimension 1, analogous to 

perimeter f uj \ 



tan"^ i-J-^i 



in dimension 2 (from j25|, (7)]) and 

, , surface 



. , (5.20) 

27T^/UJ L \kJ 2i ^ ' 

in dimension 3 (from equation in abstract]). These formulas are strikingly dimension- 
dependent, yet it must be possible to obtain them all as inverse Laplace transforms of 
the heat trace ()3.21j) . (Indeed, ()5.17j) - ()5.2()j) can be verified with the aid of the Laplace- 
transform table in 1].) The culprit, obviously, is the d-dependent power of t in (I3.2ip . When 
(|3.2ip is expanded as a power series (see (j3.22j) ). as is traditionally done, the inverse Laplace 
transforms become elementary and yield series for pK,av in powers of uo^^ whose coefficients 
depend on c/ in a relatively simple way. 

Remark: In general, p^.av for c? -|- 2 is essentially the antiderivative of p^.av for rf, as 
demonstrated by ()5.20|) and ()5.18|) . This fact can be seen either as a result of multiplying 
the Laplace transform by or as a result of convolving the eigenvalue density, with respect 
to tu^, with that of Euclidean R^, which is constant. 

To treat P and B we need the integral formula jisi (3.944.5,6)] 

^ r^-ie-""cos(o;r + 5)cir = (j-l)!(o;2 + r2)-^/^cos(j0 + 5) (^0 = tan^^ ^) . (5.21) 
We find that the inverse cosine transform of P is 

P.,per = — E ( ) (-S/^y + r'V'^ cos{2nLuJ + J0) (5.22) 

n=l j=l ^-^^ 
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and that of B is 

2 oo n+1 

P.,bdry = - y] y] c{n,j - l)/€^-i(cu2 + r^)-^/^ cos(2nLa; + (5.23) 

n=l i=2 

These formulas display p^.per and p^.bdry explicitly as sums over the periods InL of the 
periodic orbits. Finally, the complete representation of the density of eigenvalues is 

Kjbdry • 

(5.24) 

Remarks: (1) Using ()3.16|) and trigonometric identities, all the terms in principle can be 
written as cos(2r;,Lu;) or sin(2nLu;) times rational functions of uj and k. (2) The series ()5.22p 
and (j5.23p are classically divergent. To what extent one should worry about this will be 
discussed in due course. 



F. The local spectral density 

An alternative approach that has some advantages is to take the inverse cosine transform 
of (j5.7j) (with y = x) before integrating over x. The intermediate result is a local spectral 
density. The density resulting from (j5.7ap and ()5.7b|l is constant in x (namely^ (PN+PK,per)), 
so nothing new happens there. The delta-function parts of ()5.7cj) and ()5.7d|l yield a local 
spectral density 



1 

TT 



J^(-l)"+^ cos(2a;(x - (n + 1)L)) + ^(-1)"+^ cos(2cu(x + (n - 1)L)) 

^ oo 

- J2 (-l)"cos(2tu(x + nL)) (5.25) 



n=0 n=l 



TT 

n=— oo 



that is odd under x ^ L — x, so that, as anticipated, the trace of this contribution vanishes 
by virtue of cancellation of the left-hand Neumann and right-hand Dirichlet effects (cf. 
Sec. IV B|) . The rest of the terms in (j5.7c|) and ()5.7d|) give the density 



jEE(-i)"-'(")(-2M'( 

n=l 7 = 1 ^ 

[cos([2(n + 1)L - 2x]uj + j(f)) + cos([2(n - 1)L + 2x]uj + j(f))] , (5.26) 



IT 

n=l j=l 

X 

whose trace over x is 

oo n 



n=l j=l 

X [sin(2(?7, + l)Luj + — sin(2nLt<; + j 



^ oo n ^ \ 

+ ^EE(-ir' ; 

n=l j=l ' 



LU 



X [sm[2nLuj + - sin(2(n - 1) Lu + j (f))] , (5.27) 
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which "obviously" (but see Appendix [X] and Sees. IV GH V H|) simphfies to 



Pftjbou 



n=l j=l 



X 



[sin(2(n + l)Lcj + j0) - sin(2(n - l)Lcu + j0)] . (5.28) 



Presumably pK,bou = Pk,iw + Ps.bdry , though that is not obvious from the formulas (see 
Sec. EH). 



G. Alternative formulas 



Still another periodic-orbit representation of the eigenvalue density is found in Ap- 
pendix El by working backwards from the transcendental equation determining the eigenval- 
ues. One has 

p + ^S{uj) = ppois = PPois,per + PPois.bou (5.29) 

where 

L °° 
[uj) = - 1 + 2^cos(2n(Lu; + 



PPoi 



ois,per 



TT 



n=l 



1 



K 



PPois.bouii^J — — ;t 

TT tU"^ + 



n=l 



Recall from that 



K 



,.2 



sm cos 



UJ 



1 + 2^cos(2ri(Lcj + 
sin(20) 



2iu 



(5.30) 

(5.31) 

(5.32) 



so the prefactor in (j5.3ip could be combined with the phase shifts. Also, the delta function 
could be artificially written as a sum over the periodic orbits similar to ()A3|) . Forgoing those 
possibilities, however, one observes that the sums in ()5.30|) and ()5.31|) are identical except 
for a factor that depends only weakly on uj. Thus 



PPois,per(t^) = a(t^)PPois(t^), 1 - a{u) 

PPois,bou(t^) = (l - a{uj))ppois{uj), o,{uj) 



k/L 



(5.33) 



Because of its L dependence, one expects ppois.per to be the part of p coming literally from 
the periodic orbits (as opposed to the closed orbits that reverse themselves and bounce off 
a nearby boundary), hence that 



Pp- 



ois,per 



Pn + Pk. 



per ; 



(5.34) 



PPois.bou + \ ^iS^^ — PK,bou — Pfi.av + Pre.bdry • 

On the other hand, we should have 



(5.35) 



P/t,bdry PPois,bdry 



— cos(2n(Lu; + 0)) . 



(5.36) 
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The total, correct eigenvalue density p does not contain a delta function at = 0. 
Therefore, when p is computed from the formulas of Sec. IV El or from those of this section, 
a term must emerge from the other terms. We previously observed that when 

K = the compensating delta comes from the other term of Pav (|5.17p . properly interpreted. 
(Of course, in that case one can use Sec. IV Bl and never introduce the more complicated 
expressions in the first place.) When (and only when) k 7^ 0, the compensating delta must be 
hidden in the trigonometric sums, p^^per + PkM^j ^^^t approach and ppois.per + PPois,bdry 

in the other. As Appendix 1X1 shows, the presence of such a term is a delicate question, since 
it can depend on the order of the terms in the series. In contrast, the formulas of Sec. IV Fl 
do not contain any delta functions, so one must not expect Ps.pcr + Pfc.bou to contribute any 
deltas in compensation. 

However, p^.bou remains somewhat ambiguous, because of the term-ordering issue. The 
first term in (j5.27p is the contribution from those pulses (j5.7c|) that struck the right boundary, 
and the second term is from the pulses ()5.7d|) that struck the left boundary first. It is natural 
to define such sums by accumulating the terms in order of increasing path length (the same 
as the coefficient of cj, or frequency of the spectral oscillations). In (j5.26p . however, that 
ordering depends on x. There are three fairly reasonable things one could do: 

1. Simplify ()5.28|) as it stands: 

P:::^ru = -EE(-l)'"T^^(^' + ^')"'^'^i^(2Lc.)cos(2nLc.^ (5.37) 

Recall that (jOSll pairs paths with the same number of bounces from the Robin bound- 
ary: pulse 4 with pulse 7, 8 with 11, etc. 

2. Shift the index in the second term of ()5.27p by one unit, and simplify: 

pitr = ^ E L- - 1 + + + ' 

n=0 j=l ^ 

(5.38) 

Hereby paths with the same average length (as x varies) are paired: 4 with 3, 8 with 
7, etc. This choice seems to us to have the greatest physical justification. 

3. In ()5.28j) shift an index by two units relative to the other, and simplify: 

ra=0 j=l ^ 

(5.39) 

(with {2n — j + l)/n = when n = j — 1 = 0). Here all terms with the same frequency 
have been forcibly combined, so that the formula looks like a sum over periodic orbits 
alone. (The labels "pretrace" and "posttrace" refer to the timing of the balancing of 
path or orbit lengths. In both cases, unlike (j5.23|) . the cosine transform is performed 
before the trace.) 

Theorem 4 The spectral implications of Theorem are summarized as follows, modulo 
terms supported at u = in the limit: 
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(a) The local spectral density consists of the constant terms, ^5.4\ ) plus \5.2^] divided by L, 

arising from the periodic- orbit pulses (|5. 7a\ ) -^5. 7t^ , plus the terms \5.25\} and i5.2(^} 
arising from the bounce-orbit pulses \5. 7d^ - \5. 7d!^ . 

(b) The contribution of the periodic orbits to the density of eigenvalues is ^5.4\ ) plus 

or alternatively \5.30\) (the Poisson formula). 

(c) The contribution of the bounce orbits to the density of eigenvalues is ( 5.i7| ) plus \5.2'J\) 

(trace-bef ore-transform) , or alternatively any of \5/Jl]) (Poisson), ( 5.c/7| ), \5/J^) . or 
\5.39) (variants of transform-bef ore-trace) . 



H. Symbolic and numerical evaluations 

We resort now to Mathematica [29], which "verifies" (j5.34|) and (j5.36p and furthermore 
reveals that (contrast ()5.35|) ) 

posttrace tr Ar\\ 

PPois,bou = P^,bou • (5-40) 

That is, for any particular n that we have tried, machine simplification shows the identity of 
the respective terms in each of these series pairs (but general n can't be handled). Also, nu- 
merical plots of the differences of respective partial sums show nothing but numerical noise. 
(The noise is usually very small (machine-precision level) but not always: The sums over j 
in such formulas as (j5.23p are numerically unstable when uj is very small. See Fig. E^b).) 
On the other hand, partial sums of the three series (j5.37|) - (j5.39p are seen to be all different. 

When plotted, all the series show, even for very small partial sums, the development of 
delta peaks at the (square roots of the) eigenvalues of the Robin problem, 

-pL=u:'y, y{L) = 0, ^(0) = «:y(0). (5.41) 
dx^ ax 

In accordance with ()5.33|) . the periodic terms alone already give the locations of the eigenval- 
ues, and the bounce contributions are needed only to normalize the delta functions correctly. 
Examples are shown in Fig. |3 For comparison, in the case L = k = 1 we computed the 
lowest-lying eigenvalues by Newton's method to be 

cui ^ 2.0288, UJ2 ^ 4.9132, tus ^ 7.9787, cu^ ^ 11.0855, (5.42) 

and the large eigenvalues by perturbation theory to be 

" - + (^^ - 2(n -1)3.3 + (^^ZT)^ - (5-43) 

with which ()5.42|1 already overlaps well. Fig. |3] agrees well with these values. 

We now turn to the issue of what happens for u near (Fig. The empirical evidence 
from the plots is that p^.per and pK,bdry , and equivalently ppois,per and ppois.bdry , are behaving 
exactly as expected: The right-hand half of a delta function is building up at the origin to 
cancel the (correct, but misleading) —^6{uj) in the formulas of Sees. IV El and IVTjI Because 
of (j5.40p . we get the same behavior from p^^por and p^"^]^^'^'^^ ; that is, the forcible reinter- 
pretation of contributions to the local spectral density from closed but nonperiodic orbits 
as if they were contributions from periodic orbits has, by the mathematical phenomenon of 
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FIG. 4: (a) Eigenvalue density on interval 1 < w < 12 computed from 1)5. 29() . partial sum rimax = 20. 
Other formulas [ ()5.24() and analogs built from (|5.37|) . ()5.38() . and (|5.39|) ] give almost indistinguish- 
able results, (b) Same as (a) but from (|5.3U|) . the contribution of genuinely periodic orbits only. 
The plots are misleading because differences in peak height are more noticeable than differences in 
peak width. Numerical integration confirms that the peaks in (a) have almost equal, unit strength, 
while those in (b) do not, in accordance with ()B5|) and ()5.33|) . (c) Eigenvalue density on interval 
50 < u; < 60 computed from 1)5. 29p with rimax = 4. (d) Same as (c) but with nmax = 20. In all 
plots, L = K = 1. 



Appendix El introduced a truly spurious term +\5{uj) into the formulas of Sec. IV Fl On 
the other hand, both p^f^^^T'^ ^'^d pJJ'bou converging to the right answer for the context of 
Sec. IV Fl Their partial sums differ, but the difference (apparently) goes away in the limit. 
These series are yielding negative delta functions that precisely cancel the delta behavior of 
PK,per , so that no 5{ijj) is ever visible in the total density (see Figs. ini^c,d)). 



I. Broader implications 

This simple one- dimensional model is trivial in comparison with many problems treated 
by similar methods in the contemporary literature. It cannot manifest quantum chaos; 
it doesn't require stationary-phase approximations; it isn't even relevant to the sort of 
semiclassical approximation that becomes necessary when a potential is added to the one- 
dimensional Schrodinger equation ()5.4H) [l^. Nevertheless, it is instructive. The mere fact 
that so much information can be obtained exactly means that the model can be understood in 
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FIG. 5: (a) Eigenvalue density on interval < u; < 3 computed from (|5.29|) (the Poisson formula) 
with riinax = 20. Cancellation of the delta function at a; = is clearly visible, (b) Difference 
between /9K,bdry (|5.23() and ppois.bdry (|5.36|) on the interval < tj < 4, rimax = 6, showing only 
roundoff error, (c) Same as (a) but computed from the pathlength-balanced formula 1)5. 38() [plus 
(|5.4|) and (|5.22|) ]. The delta function is now missing, (d) Same as (c) but computed from the "naive" 
formula (|5.37j) : the difference presumably is purely truncation error. In all plots, L = k = \. 



complete detail and stands as a benchmark against which partial and approximate solutions 
of more complicated models can be appraised. 

In particular, the delicate problem of the delta function at = is likely to have broader 
implications. In more complicated problems such precise analysis at the bottom end of the 
spectrum is usually lost in the noise of the semiclassical and stationary-phase approxima- 
tions. Periodic-orbit reconstructions of spectra are observed to be surprisingly good at low 
frequencies, except very close to 0, where a spurious peak often occurs. (Look at Fig. 17 
of 1^ and Fig. 29 of |l3|.) Our analysis suggests that this phenomenon represents not a 
breakdown of the semiclassical approximation so much as an ambiguity in the ordering of 
the terms in the badly convergent periodic-orbit sum. 

A related point is the significance of contributions from orbits that are closed but not 
periodic. Usually such contributions are discarded in the process of stationary-phase ap- 
proximation. However, they are needed to get the boundary terms in pav fsl, la], and Jaffe 
and Scardicchio (20|, |2J] have recently emphasized their importance in calculations of total 
vacuum energy. As semiclassical calculations are carried beyond the lowest order in ^, these 
orbits will need to be included, and one must grapple with the question of whether they must 
be kept separate, as in (j5.27p . or will be effectively absorbed into the periodic contributions 
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as in pK,bdry and p^'bou^'^*^ In our model they did not influence the location 

of the eigenvalues, but they were needed for proper normalization; is there some broader 
significance to that? 
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APPENDIX A: NONABSOLUTE CONVERGENCE OF OSCILLATORY DIS- 
TRIBUTIONS 

In the periodic-orbit approach to spectral theory, sums of the type 

oo M 



am(n, uo) cos(nco' + 5m) (Al) 



n=0 m=l 



are often encountered. Typically M is a fixed, small number (such as 2), and am{n, uo) varies 
slowly with uj but may become singular as — and, most important, does not approach 
as n — *• oo. Such a series does not converge at all according to the definitions of classical 
analysis, but nevertheless it may converge in the sense of distributions. The question then 
arises whether the distributional limit can depend on the order in which the terms are added, 
in analogy with numerical series that are only conditionally convergent. We consider here 
distributions defined on the nonnegative real line. 

Theorem 5 In \A1\) assume that M is independent ofn and that the functions am{n, u) are 
smooth (C°°) for < < oo and they and all their derivatives are uniformly polynomially 
bounded as n ^ oo. Then any two orderings of the terms in (Al) define distributions that 
coincide on test functions with support in the interior (i.e., 4>{u!) = in a neighborhood of the 
origin). As distributions on test functions defined in [0, oo) they may differ by a distribution 
supported at u = (necessarily a linear combination of 5{uj) and its derivatives). 

Proof: The definition of distributional convergence is that for any test function 0, 

oo M „oo 

/ am{n, u) cos{nuo + 5m)4'{'^) duj (A2) 

n 1 Jo 



n=0 m=l 



converges in the classical sense and defines a continuous linear functional of 0. If (f){uj) = 
in a neighborhood of the origin, then repeated integration by parts. 



Jo 



oo 



1 d^P 

am{n,uj)cos{nuJ + Sm)(l){uj)duj = / — -^-^ cos{nuj + Sm)am{n,uj)(j){uj) du 



J cos{nuj + 6m)-^^[amin,uj)(f){uj)]duj 
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shows that each integral in ()A2|1 falls off faster than any power of n. Taking p sufficiently 
large guarantees that ()A2|) converges absolutely and hence can be reordered at will. If 
the support of includes the origin, the endpoint terms from the partial integrations may 
converge to a well-defined delta-type distribution. If not, it may still be possible to extend 
the distributions to such test functions by "regularization" ; in that case the coefficients 
in the delta sum are somewhat ambiguous. 
Example: Consider 



— [sin((n + 1)0;) — sin(nco')]. (A3) 

n=0 ^ 

From one point of view, shifting the index in the first term yields 

— [sin(nco') — sm(nuj)] = 0. 

n=0 ^ 

On the other hand, the nth partial sum of ()A3|) as written is 

N-l ^ ^ 

E — [sin((r;, + 1)^^;) — sm(nuj)] = — sin(A^co'), 
r,i I,} 



n=0 



and the latter converges to ^6{uj 



2 

sm{Nuj)(l){uj) dcu = / - sm{Niu)[(f){0) + 0{uj)] duj 



UJ Jo ^ 

coo 



= 0(0) / -sm{z)dz+ [ sm{Nuj)0{l)duj; 
Jo z Jo 

the second term vanishes as iV ^ 00 by the Riemann-Lebesgue theorem, and the integral 
in the first term equals f [11, (3.721.1)]. 

Remark: That the series rearrangements in the body of the paper obey the polynomial 
boundedness requirement in Theorem El follows from the bound ^j, (2.14.13)] 

\L\_,{x)\e-''/^ <n (x>0) (A4) 

on the Laguerre polynomial. Individual terms in ()5.28|) (for instance) with J ~ f can grow 
exponentially with n, but we have never realigned terms with differing values of j. 



APPENDIX B: RECONSTRUCTION OF THE PERIODIC-ORBIT SUM BY 
POISSON SUMMATION 

The following is a close analogue of a two-dimensional calculation in Sec. 2.2 of jl^. 
It requires knowledge of the eigenvalue condition (jBlll (but not of explicit formulas for its 
solutions). 

The square roots of the eigenvalues of the Robin-Dirichlet problem ()5.4H1 with k, > are 
the positive roots of 

tan(La;) = - - , (Bl) 

K 
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which may be parametrized as 

fc7r = Ltu + 0, = tan-^(^^j, k=l,2,.... (B2) 

One can view ()B2|) as giving as a function of u, which naturally extends to the whole real 
line as a monotonic and odd function with derivative 

aUJ TT \ CU"^ + / 

Thus Uk is defined for all integers, with u_k = ujk and uq = (which is not an eigenvalue). 
So the eigenvalue density is 



oo _ oo _ 

J2S{uj-u;,) = - J2 S{u;-\u,\)--6{uj). (B4) 

k=l k=—oci 

By the Poisson summation formula il9] (and ()B3|) and ()B2|) ) one has 
- ^ 5(0; - = - ^ / rffc5(a; - \uj,\)e'-^'- 

k=-oo n=-oo''~°° 

n=-oo ^ 

°° 1 / \ 

= 0(^) -f^ + )cos(2n(L.; + </))), 

n=— oo ^ ' 

since = |ci)| has two roots u for u > and none for u < 0. Therefore, 

oo ^ 0(uj) oo ^ K \ 

y 5(cu -ujk) = - -S{uj) + L + + 2 y ( L + ) cos(2n(Ltu + 

k=l L n=l ^ ^ 

(B5) 

It is natural to associate the terms in ()B5|) with prefactor L directly with the periodic 
orbits and to regard the other terms as the traces of the contributions of the bounce orbits. 
With this interpretation, and confining attention to > 0, we have arrived at ()5.29p - ()5.3ip . 
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